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On  Calculating  with  Lower  Order 
Chebyshev  Splines 

Mladen  Rogina  and  Tina  Bosner 


Abstract.  We  develop  a technique  to  calculate  with  Chebyshev  Splines 
of  orders  3 and  4,  based  on  the  known  derivative  formula  for  Chebyshev 
splines  and  an  Oslo  type  algorithm.  We  assume  that  splines  in  the  reduced 
system  are  simple  enough  to  calculate.  Local  bases  of  Chebyshev  splines  of 
order  3 and  4 can  thus  be  evaluated  as  positive  linear  combinations  of  less 
smooth  Chebyshev  B-splines.  The  coefficients  in  such  linear  combinations 
are  discrete  Chebyshev  splines,  normalized  so  as  to  make  a partition  of 
unity.  There  are  a number  of  interesting  special  cases,  such  as  Foley’s 
i/-splines,  Chebyshev  polynomial  splines  (g-splines),  and  splines  in  tension 
which  can  be  calculated  stably  by  such  formulae. 


§1.  Introduction  and  Preliminaries 

It  is  an  important  fact  in  the  univariate  polynomial  spline  theory  that  splines 
can  be  represented  as  linear  combinations  of  compactly  supported  basis  func- 
tions, which  we  can  calculate  in  various  ways  by  stable  and  fast  numerical  algo- 
rithms. This  can  be  extended  to  some  other  well-known  spaces  of  splines,  such 
as  trigonometric  and  hyperbolic  splines,  where  nice  three-term  recurrences  of 
de  Boor-Cox  type  exist;  this  also  applies  to  a less  interesting  case  of  Chebyshev 
splines  with  equal  weights.  These  issues  have  been  discussed  in  [12],  where  we 
can  also  find  a negative  result  concerning  the  existence  of  such  three-term  re- 
currence relations  in  general.  Polynomials,  however,  form  an  algebra,  while  in 
other  cases  multiplicative  properties  are  replaced  by  other  algebraic  formulae, 
such  as  addition  formulae  for  trigonometric/hyperbolic  functions  or  similar 
identities. 

We  shall  say  that  an  interval  [a,  6]  is  measurable  with  respect  to  the 
measure  vector  da  :=  (da2;  • . . dan)T  if  it  is  measurable  with  respect  to  the 
positive  Lebesgue-Stieltjes  measures  dcr,,  i — 2, . . .n.  Then  for  x € [o,  6]  we 
can  define  generalized  powers  (Chebyshev  system)  (1,  u2, . . . un}: 

j>X  rX  fin- 1 

u2{x)  = / d<T2(t2);  ■ ■ • ; un{x)  = / dcr2(t2) ...  dcrn(tn).  (1) 

J a J a J a 
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If  all  of  the  measures  dai  are  dominated  by  the  Lebesgue  measure,  then  they 
possess  densities  1-,  i = 2, . . . n;  if  pi  are  smooth,  ^ € Cn_,+1, 

then  {1,«2, ...«n}  is  called  an  Extended  Complete  Chebyshev  System  (ECC- 
system),  referred  to  as  ECT  system  in  [11]).  Further,  we  shall  assume  that  such 
an  integral  representation  has  been  found,  and  the  measures  for  the  Chebyshev 
space  determined  in  such  a way  that  we  know  the  generalized  powers  explicitly. 
This  may  not  always  be  easy,  and  the  choice  of  measures  may  not  be  unique. 
It  is  known  that  Chebyshev  B-splines  that  make  a partition  of  unity  exist 
in  this  general  case,  and  even  that  abstract  recurrences  [5,1]  resembling  the 
polynomial  ones  exist.  Almost  nothing  can  be  said  about  their  numerical 
stability,  at  least  not  until  we  employ  a special  measure  vector,  whence  the 
abstract  construction  gets  difficult. 

Recently,  other  techniques  based  on  blossoming  have  been  found  for  the 
Chebyshev  splines  [6],  which  are  more  promising  so  far  as  evaluation  and 
numerical  stability  is  involved.  In  case  of  polynomial  weights,  one  obtains 
Chebyshev  polynomial  splines  [7],  though  here  again  the  underlying  algebraic 
properties  of  polynomials  are  implicitly  used. 

In  Section  2 we  give  some  formulae  for  Chebyshev  systems  of  orders  3 and 
4 which  express  locally  supported  splines  (being  piecewise  in  these  spaces) 
as  positive  linear  combinations  of  less  smooth  splines  in  the  same  space.  We 
will  see  that  the  coefficients  are  related  to  integrals  of  splines  in  the  reduced 
system,  which  is  defined  to  be  a Chebyshev  system  like  (1)  corresponding  to 
the  reduced  measure  vector,  that  is,  for  each  6 C [a,  6] 

d<r^(6)  :=  (dai+2(S),. . . ,dan{6))T  € Rn~{i+1),  t = l,...n-2. 

If  S(n,d<r)  :=  span{l,  «2, . . . «n},  then  the  generalized  derivatives  Lj^cr  • — 
Dj  ■ ■ ■ D\Dq,  where 


Djf(x)  :=  lim 


f(x  + 6)~  f{x) 
s^o+  aj+i(x  + 6)  - 0j+ i(x) 


j = 


(2) 


are  linear  mappings  S[n,dxr)  — > S(n  — j,da^).  For  a partition  A = of 

an  interval  [a,  6],  and  a given  multiplicity  vector  m = (nj, . . . ,rik)T , (0  < ri{  < 
n),  we  shall  denote  by  (fi . . . t2„+k}  [11]  an  extended  partition  in  the  usual  way 
(see  Fig.  1 and  Fig.  2): 


t\  = • • • = tn  = x0, 

tn+fc+l  = • ' • = t2n+k  = xk+ 1) 

^ ^ ^n+k  ~ Xk , • • • , Xk  - 

' V " ' V ' 

ni  nk 

S(n,m,d(r,  A)  is  the  spline  space  spanned  by  functions  being  piecewise  in 
S(n,da)  [11].  Chebyshev  B-splines  in  S(n,m,dcr,  A)  have  compact  support 
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[t,,ti+n\,  and  we  shall  henceforth  assume  that  these  are  unique  such  splines 
such  that 

n+K 

(*)  = !;  i = 0, . . . n — 2,  (3) 

•=i 

where  K :=  Ylini-  We  aim  to  show  that  for  lower  order  Chebyshev  splines, 
one  can  express  Chebyshev  B-splines  as  linear  combinations  with  positive 
coefficients  of  B-splines  with  multiple  knots.  Hopefully,  these  can  be  computed 
efficiently  by  some  interpolating  formula  or  otherwise.  In  Section  3 we  give 
some  examples  of  how  the  theory  can  be  used  to  construct  Chebyshev  B- 
splines  in  some  known  polynomial  cases,  like  weighted  splines  and  g-splines, 
but  also  in  the  non-trivial  case  of  some  other  useful  Chebyshev  splines.  To 
achieve  this,  we  need  some  technical  results.  One  is  the  derivative  formula, 
stating  that  for  x € [a,  6],  and  a multiplicity  vector  m whose  components 
satisfy  n*  < n — 1 (i  = 1, . . . k),  the  derivative  of  a T-spline  is 


Li,dtrT^d<T(x)  = 


Cn-lii) 


t;TiU(*) 
Cn-l(i  + l)' 


(4) 


where 

pti+n  — l 

Cn-i(i)--=  l 2.  (5) 

For  continuous  a,  the  proof  is  similar  to  that  in  [9];  a somewhat  longer  proof 
involving  only  determinant  identities  exists,  and  relies  only  on  the  fundamental 
theorem  of  integral  calculus  [11]: 


f(b ) - /(a)  = f Lhdaf{t)do2(t), 
J a 


(6) 


which  holds  under  very  weak  hypothesis  on  the  measures. 

In  certain  cases  one  may  construct  splines  defined  on  triplets  of  knots 
(Lagrangian  splines)  by  a generalized  Taylor  formula.  We  give  a variant  of  the 
Taylor  formula,  amenable  to  generalization,  that  can  be  useful  for  lower  order 
Chebyshev  systems: 

Lemma  1.  If  f and  da  = (da2,  ■ ■ ■ , dos)7  are  such  that  L*/  : = Liydaf  exist 
and  are  measurable  with  respect  to  dcri+l!(i  = 1, ...  4 ) on  [a,  5],  then 


f(x)=f(a)  + Lif(a)[  dcr2(s2)  + L2f(a)  f da2(s2)  f da3{s3) 
J a J a J a 

rx  pS2  ps 3 

+ L3f(a)  / da2(s2)  / da3{s3)  / d<74(s4) 

J a J a J a 

px  ps  2 pS3  ps  4 

+ / dcr2(s2)  / da3(s3)  / da4(s4)  / L4f(s5)da5{s5). 

J a J a J a J a 
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Proof:  The  hypothesis  and  (6)  enable  us  to  use  standard  arguments  of  math- 
ematical analysis  in  developing  a Taylor  series  expansion 


f(x)  =f(a)+  j (Lif(s2)  - Lif(a)  + Lif{a))da2{s2) 
J a 

-f(a)  + Lif(a)  j da2{s2)  + etc.  □ 

J a 


Lemma  1 is  in  fact  a way  of  writing  a Taylor  expansion  for  L-splines, 
see  [13],  pp.  425-426. 


§2.  A Knot  Insertion  Algorithm 

We  shall  henceforth  assume  that  Chebyshev  B-splines  T.3d_(1)  can  be  evaluated 
at  the  knots  in  a numericaly  stable  way.  This  is  a sound  hypothesis,  since 
by  (4), 


T''da{x)  ~~mL  T'’da(i)da2  ~ WTT)  Jti+l 

If  we  suppose  that  the  multiplicity  vector  m is  such  that  T?da  6 C\a , 6]  (that 
is,  rii  < 2),  then  we  have 


and 


J ft  >+2  J fti+3 

Tlda{ti+2)  = 1-C2(i  + 1)  /(  +i  Th,da^d<J2  = Cz(.  + 1)  /fj+a  Ti+i 


Since  the  construction  of  two-interval  supported  “linear”  splines  is  easy,  sta- 
ble evaluation  of  T?da(ti+i),  Tfda.(U+ 2)  amounts  to  finding  an  integration 
formula,  preferably  of  Gaussian  type.  The  important  thing  is  that  in  this 
case  we  do  not  have  any  dangerous  subtractions  potentially  leading  to  large 
floating  point  errors,  as  in  (4). 

Theorem  2.  Let  T3dff(1)  g 5(3,  m,  da^,  A)  be  a Chebyshev  3rd  order  B- 
spline  associated  with  the  multiplicity  vector  m = (1, . . . 1)T,  and  let  us  as- 
sume that  T3dff(I)  € 5(3,m.  da1'11),  A)  are  Chebyshev  B-splines  associated 
with  the  multiplicity  vector  to  = (2, . . . 2)T  (Fig.  1)  on  the  same  partition.  If 
{fi, . . . tfc +6 } and  {f i , . . . t2fc+e}  are  the  associated  extended  partitions,  and  r 
an  index  such  that  tj  = t r < tr+i,  then  for  j = 1, . . . k + 3, 


= Tj  d<r<.i)(tj+i)TT  d 
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Fig.  1.  Double  knots. 


Proof:  Let  Tf  T3rfff(1) . Since  Tj(tj)  = T?(tJ+3)  = 0 and  the  same 
holds  for  the  first  generalized  derivatives,  we  conclude  that  two  out  of  three 
coefficients  representing  T3  on  each  interval  of  its  support  are  zero.  For 
x € (tj,tj+i)  we  have  that  Tf(tj+ 1)  = 63(r)T'3(£j+i).  Since  (3)  applied  to 
<S(3,m,  dcr(l\  A)  implies  that  T3(tj+i)  = 1,  it  remains  to  show  that  the  mid- 
dle coefficient  equals  1.  For  x € (tj,tJ+i),  property  (3)  again  leads  to 

Tj-i(x)  + Tj(x)  + Tj+i(x)  = 1. 

We  expand  T3,  i € {j  — 1 ,j,j  + 1}  in  <S(3,m,  da^,  A),  and  rearrange  the 
terms  to  obtain 

1)  + Tf(tj+ 1)]  + fi?(r  + l)T3+1(x) 

+ ^H-2(x)[r/(fi+2)  + Tf+2{tj+2)\  =1. 

The  expressions  in  square  brackets  are  equal  to  1 by  (3)  applied  at  the  knots 
tj+i,tj+2-  But  the  partition  of  unity  property  must  also  hold  for  Chebyshev 
B-splines  in  <S(3,m,  dall\  A),  and,  since  the  expansion  of  unity  in  this  space 
must  be  unique,  we  have  Sj(r  + 1)  = 1.  □ 

The  more  important  “cubic”  version  follows  from  Theorem  2 and  the 
derivative  formula  (4): 

Theorem  3.  Let  Tjda  £ <S(4,m,  da,  A),  Tyda  £ 5(4,  m,  da,  A),  and  let  m,  m 
be  multiplicity  vectors  as  in  Theorem  2.  Then  there  exist  positive  fif(i),  de- 
pendent on  da 2,  such  that  T^da  = YH=r  where  r = rj  satisfies  tj  = 

tTj  < tP  +1.  If  the  the  extended  partition  is  {fi, . . . tfc+s},  and  {ti, . . . t2fc+s} 
is  the  extended  partition  with  double  interior  knots,  then  6j(i),  i = r, . . . r + 3 
are  explicitly  determined  by 


s,()  = 

j Tj,da(  + &(r  + !)  + ^ydo-(o(it+2)C'(r  + 2) 

rV)(*j+i)C(r)  + C(r  + l) 

+ ^ + 1)  + ^,,(<,+2)^  + 2)’ 

T*+hdaM(tj+3)C(r  + 4)  + C(r  + 3) 

+ 2?+1.a»o)(^+2)^(»'  + 2)  + C-(r  + 3)  + + 4)’ 


Sj(r  + 3) 


Tj+i,dgO)(tj+3)C{r  + 4) 

T?+  moo.  (tj+2)C(r  + 2)  + C(r  + 3)  + T3^^,  (tj+3)C(r  + 4)  ’ 
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where  as  in  (5), 

£(*')=  / Tlda(r,d02.  (7) 

support 

Proof:  We  expand  Tjda  in  terms  of  less  smooth  T*da: 


r+6 


r-t-o 

TfM*)  = E «/(0£ 


•fxh 


J/T\X 


i=r—2 


and  apply  (4)  to  obtain 

Tf'do{x)  T?+l'd(r(x)  _ Sj(i)  - 6j(i  - l)-3 

c3(j ) CsO’  + i)  V <53(0 

We  may  then  use  Theorem  2 to  expand  Tjd(T,  and  write  a linear  system 
for  6j,  which  has  the  above  explicit  solution.  Details  are  omitted  since  the 
construction  very  much  resembles  the  construction  of  (/-splines  in  [10].  □ 

It  is  not  difficult  to  prove  that  6j  are  discrete  Chebyshev  splines,  that  is, 
they  form  a partition  of  unity,  i.e.,  t>iU)  — 1 [10];  this  also  holds  for  general 
order  splines. 

We  also  note  that  coalescence  of  the  knots  yields  an  expression  for  a 
complete  “Chebyshev  cubic”  spline  on  triple  partitions. 

§3.  Examples 

We  investigate  how  much  of  the  above  theory  can  be  applied  to  some  known 
spaces  of  Chebyshev  splines,  and  what  special  properties  must  be  used  to 
obtain  the  stability  of  the  algorithm  for  the  evaluation  of  Chebyshev  B-splines. 

Remark  4.  Minimizing  elastic  energy  of  an  inhomogenuous  rod  leads  to  a 
minimization  problem  for  the  functional 

(E(s)I(s)w,,(s))2  ds  — *■  min, 

where  E is  Young’s  modulus  of  elasticity,  and  I is  the  moment  of  inertia  of  its 
cross  section.  We  think  of  the  rod  as  being  made  of  pieces  of  different  material 
(E  piecewise  constant),  or  different  cross  section  (I  is  piecewise  constant)  on 
intervals  [x,,Xi+x)  that  partition  [a, 6).  In  either  case,  the  Euler  equation  is 
I/4,dcrw  :=  ( wu ")"  = 0,  where  da  :=  (dx,da,dx)T),  and  da  is  the  measure 
generated  by  the  piecewise  constant  density  w : = 1 /El,  = : w, . 

In  CAGD  such  splines  are  known  as  splines,  introduced  by  Foley  [2]. 
It  follows  from  the  variational  formulation  that  (/-splines  possess  continuous 
derivatives,  and  that  jumps  in  second  derivatives  must  be  continuous: 


Wi+iu”(Xi)  — Wiu"(Xi)  = 0. 


(8) 
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If  the  breakpoints  for  w are  points  of  the  partition  A,  it  is  not  difficult  to  see  (2) 
that  (8)  is  equivalent  to  the  continuity  of  the  second  generalized  derivatives 
1,2, da  across  the  knots.  We  may  evaluate  Tfda  by  scalar  products  of  positive, 
known  quantities;  for  instance  it  follows  (with  the  knots  enumerated  as  in 
Theorem  3),  that 


rp4  72  (4)  04  , 7l(4)  +72(5)  n4  . 71(7)  +73(8)  4 . 7|(8)  f, 4 

J-0  — n -a. i ■d4  “I  ii  3ii  "I  Ti~^n  n6  "r  ■■  ®“  -D 


1172  II 


hall 


hill 


hill 


75 


where 

72 (4)  = - fa),  72(5)  = U ~ t3, 

72(7)=t5-t4,  72 (8)  = dg(^'’'tg)(*6  ~t4)’ 

hill  = 1^7i3(*)  for  1 = 2,3. 

i 

One  should  note  that  Bf  are  ordinary  polynomial  splines  (by  raising  the  mul- 
tiplicity of  the  knots  we  avoid  the  second  derivative  condition),  and  thus  are 
readily  calculated  by  the  de  Boor-Cox  recurrence.  This  is  an  example  of  a 
Chebyshev  system  which  is  not  an  ECC-system,  since  da  is  generated  by  a 
non-continuous  density. 

Remark  5.  If  da  in  Remark  4 is  a Lebesgue  measure,  the  standard  knot 
insertion  formula  for  a “homogenuous  rod”  i.e.,  cubic  splines  appears: 


4 _h  - f2  jj4  + *4  - h £4  fg  -U  £4  t6  -f4  -4 


ts  ~ h 


h — f2 


BZ  + 


- h 


— h 


-By. 


Remark  6.  In  the  last  polynomial  example,  we  consider  the  g-splines  intro- 
duced by  Kulkarni  and  Laurent  [3] , which  are  piecewisely  linear  combinations 
of  the  functions  in  the  canonical  system  {l,u2,u3,u4}: 


dt  2, 


u2{x)  = f 

J a 

fX  rt2 

u3{x)=  / dt2  / q{t3)dt3, 

J a J a 

rx  /*<2  /»<3 

u4(x)  = / dt2  / q(t3)dt3  / dt4, 
J a J a J a 


where  q is  piecewise  linear: 


q(x)\lti,u+1y.=  ^^(x-U)  + qi, 


hi  : = ti+ 1 — U,  and  Qi  > 0.  This  may  be  thought  of  as  an  elastic  rod  with 
elastic  properties  changing  in  a reciprocally  linear  way;  the  more  physically 
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Fig.  2.  Triple  knots. 


sound  model  in  which  material  properties  are  changing  linearly  does  not  lead 
to  “Chebyshev  polynomial  splines”  [7],  but  involves  logarithmic  weights,  and 
does  not  seem  to  be  very  useful. 

Some  parts  of  the  construction  can  be  made  more  explicit,  i.e.,  (7)  can 
be  integrated  to  obtain 


C3(r  - 1) 


+ qi  + 2qi- 


4 

/y  / \ 1 Qi^i  i u i u i 17i+2<b  + l 

Cs(r)  =-  — — - + ft,  + ftt+i  + 

4 + 2?i+i  2g;+i  + 9i+2. 


We  can  stably  calculate  “parabolic”  B-splines  required  by  Theorem  2 in  two 
important  points  via 


T?(*i+i)  = 


hj{qi+2qi+1) 

6 C2(i) 


rp3,.  \ _ hi+2{2qi+2  + ql+ 3) 

Ti  [U+2>  6 C2{i  + 1) 


where 


B?(t)q(t)dt 


e[(?:  + 2<7i+i)hi  + {2qi+i  + <?j+2)hj+i]. 

b 


One  can  in  fact  express  T?da  in  terms  of  B-splines  on  a triple  net  (Bernstein 
polynomials).  The  enumeration  of  the  knots  is  as  in  Fig.  2,  with  quadruple 
knots  at  both  ends  and  triple  knots  ti  in  the  interior.  If  s is  an  index  such 
that  ti  — tg 2 — tg — r — tg  ^ , then 


ji4  _ 1 2gj  hi  -5  1 / hi  T <b+i 

1 C3{r  —1)  2gj+g,+i  4 s 2 C3(r)  ' 4 

qi+2  . 1 9i+2  hi+ 1 

2gi+i+gi+2  4 / s 1 C3(r)  "^qi+i  + qi+2  4 s’ 

and  an  analogous  result  holds  for  T*. 

Remark  7.  Let  us  now  see  how  the  theory  may  be  applied  to  the  “real” 
Chebyshev  case  by  considering  tension  splines  with  uniform  tension  parameter 
that  are  piecewise  in  the  null-space  of  the  differential  operator  D2(D'2  — p2), 
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where  p > 0 is  the  tension  parameter;  see  [4]  for  some  more  recent  references. 
We  may  factor  the  differential  operator  as 

D(D  + ptanhpx)(D  — ptanhpx)D  = ( — ^ — £>)(cosh2  pxD)(— — D)D, 

cosn  px  cosn  px 

and  identify  the  measure  vector  da  = (ds2,  coshps3ds3,  Other 

measure  vectors  may  be  used,  but  this  has  the  advantage  that  L\tda  = D 
maps  the  tension  spline  space  to  the  space  of  hyperbolic  splines,  where  nice 
recurrencies  exist.  To  apply  Theorem  3,  we  need  B-splines  in  the  reduced 
system  [14]: 


x £ ( tj , fj+i), 


X £ (tj+l , fj'+2)> 

X £ (tj+ 2)  3)) 

where  the  normalisation  constant  C : = cosh  |(tJ+2  — tj+ 1)  ensures  the  parti- 
tion of  unity  (3).  The  problem  of  how  to  calculate  T*da  remains.  If  we  use 
techniques  like  in  Remark  6 to  express  tension  splines  on  triple  nets,  we  finally 
arrive  at  the  numericaly  nasty  formula  for  a T-spline  with  support  [U,U+ 1], 
(h  : = ti+ 1 - tj),  that  is  is  given  by 

(sinhph  — ph)(coshp(x  — ti)  — 1)  — (cosh  ph  — l)(sinhp(x  — t{)  — p(x  — ti)) 
2(sinhp/i  — ph)(ph/ 2 coshp/i/2  — sinhpfi/2) 

The  above  formula  is  a special  case  of  the  integrated  version  of  the  derivative 
formula  (4).  Taylor  expansions  may  be  used  to  calculate  the  above  expression 
for  a small  p (approx,  p < 0.5  in  double  precision,  IEEE  standard),  and  also 
an  asymptotic  formula  for  the  ultimate  almost  linear  spline,  ideas  similar  to 
the  ones  used  by  Rentrop  [8].  In  the  middle  range  it  is  best  to  utilize  the  gener- 
alized Taylor  expansion  from  Lemma  1 in  the  vicinity  of  the  inflexion  point  of 
Tidffi  which  is  defined  as  the  solution  of  the  equation  Lz,daTida{a inflex)  = 0. 
On  the  standard  interval  [0, 1],  one  obtains  a;nflex.  = ^ log  p!^exp’(-py-1i  i both 
limits  never  approach  boundaries. 

For  h = 1 we  can  thus  obtain  an  absolutely  stable  formula  in  the  range 
0.5  < p < 700  with  very  few  arithmetic  operations.  Translation  invariance 
of  splines  in  tension  enables  calculation  on  any  interval.  We  also  note  that 
in  order  to  have  a complete  algorithm,  one  must  have  a kind  of  Gaussian 
integration  formula  in  closed  form  to  calculate  the  normalisation  constants  (5). 
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sinh2  (p/2)(z— t,) 

3inh(p/2)(t3+2-tj)sinh(p/2)(tJ+i-tj)’ 

sinh2  (p/2)(x—tj ) sinh2  (p/2)(t,+2-i)  . 

m 3 — pJ  sinh(p/2)(L+2-L)sinhCp/2)(L+2-L  + i)i_ 

i'daW  sinh2  (p/2)(tj+3-i)sinh2  (pl2)(x-tj  + ix) 

sinh  (p/2)(tj+3-t3+i)sinh  (p/2)(t3+2-*i+1)  > 

sinh2  (p/2)(t,+3  — s) . 

. sinh (p/2)(t3+3— tj+l)sinh{p/2)(tj+3— q+2)’ 
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